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We present numerical calculations of the equation of state for dense matter in high 
magnetic fields, using a temperature dependent Thomas-Fermi theory with a magnetic 
field that takes all Landau levels into account. Free energies for atoms and matter are 
also calculated as well as profiles of the electron density as a function of distance from 
the atomic nucleus for representative values of the magnetic field strength, total matter 
density, and temperature. The Landau shell structure, which is so prominent in cold 
dense matter in high magnetic fields, is still clearly present at finite temperature as 
[ — | long as it is less than approximately one tenth of the cyclotron energy. This structure is 

reflected in an oscillatory behaviour of the equation of state and other thermodynamic 
■ properties of dense matter and hence also in profiles of the density and pressure as 

functions of depth in the surface layers of magnetic neutron stars. These oscillations 
2 ■ are completely smoothed out by thermal effects at temperatures of the order of the 

cyclotron energy or higher. 
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1. Introduction 

The temperature of a neutron star exceeds 10 11 K at its birth in a supernova core collapse. 
Initially the star cools very rapidly due to a copious emission of neutrinos (and anti-neutrinos) 
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but once the interior temperature has become of the order of 10 9 K the cooling becomes slower 
although it is still dominated by neutrino emission for thousands of years. At later stages when 
the temperature has fallen sufficiently the emission of electromagnetic radiation from the surface 



becomes the dominant cooling mechanism (see e.g. Prakash et al. 1997; Pethick 1992; Nomoto & 
Tsuruta 1987j ; |Umeda et al. 1994Q . 



An important ingredient in calculations of neutron star cooling is the thermal structure of the 
envelope, i.e. the surface layers and the outer crust of the star. Although the envelope contains 
only a very small fraction of the star's mass, it controls the heat transport from the high density 
interior to the surface and thus determines the relationship between the surface temperature and 
the temperature of the stellar core, which is isothermal for sufficiently old stars. Among other 
things this relationship plays a crucial role in relating results of theoretical cooling calculations to 
observations of radiation from neutron star surfaces (see pudmundsson et al. 1985| , and references 
therein). 

If the neutron star has a high magnetic field the properties of matter close to the surface are 
strongly modified by the field as compared to the zero field case, thus changing both density and 
temperature profiles in the envelope (see e.g. Rognvaldsson et al. 1993] (hereafter Paper I), and 
references therein). The field may also have a considerable effect on the spectrum of radiation 



emitted from the stellar surface as well as its polarization and angular distribution (see e.g. Pavlov 
et al. 1995| ). 



In order to investigate the thermal structure of neutron star envelopes one needs as physical 
input both the opacity as well as the equation of state of matter for the relevant temperatures 
and densities. The case of iron envelopes in zero magnetic field was investigated in detail by 



Gudmundsson et al. (1983) , and more recently the thermal structure of envelopes with accreted 



layers of lighter elements has been investigated by Chabrier et al. (1997) and Potekhin et al 



1997). Several authors have also included the effects of strong magnetic fields to varying degrees 



of accuracy (e.g. Hernquist 1985; Van Riper 1988; Schaaf 1990a, b; Page 1995). The paper by 
Yakovlev & Kaminker (1994) gives a comprehensive review of the physical properties of neutron 
star envelopes with magnetic fields, including a general discussion of magnetic opacities. More 



recently [Potekhin (1996)| and |Potekhin &: Yakovlev (1996)| have published practical expressions for 
electric and thermal conductivities along magnetic fields under conditions occurring in neutron 



stars, and Lai fc Salpeter (1997) have investigated the effects of different phases of hydrogen on 
the properties of the atmospheres and surface emission of strongly magnetized neutron stars. 

In this paper we extend our previous work on the zero temperature equation of state of dense 
matter in high magnetic fields (Paper I) by including the effects of temperature for all Landau 
levels. As before we use the pure Thomas-Fermi (TF) method for matter in a magnetic field, 
where by pure we mean that correction terms such as the exchange (Dirac) correction and the 
von Weizsacker gradient correction are not included in the calculations. For the problem under 
discussion this is in fact a good approximation as emphasized by [Fushiki et al. (1992) , who present 
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a detailed treatment of magnetic TF theory and give extensive references to earlier work (see also 
Paper I and Lieb et al. 1994| ). 

The paper is organized as follows. In §2 we first give a short general overview of TF theory 
at arbitrary temperature and then specialize the discussion to hot matter in high magnetic fields 
with all Landau levels taken into account. There we also discuss the range of validity of the theory 
and its scaling properties. Our numerical methods are described in §3, and in §4 we present the 
results of the calculations for the equation of state and the sound velocity of hot neutron star 
matter as well as isothermal density profiles of surface layers of neutron stars. As in Paper I we 
have chosen iron, with Z = 26 and A = 56, as our reference element in the numerical calculations, 
and the range of magnetic field strengths covered is from 10 10 to 10 13 G. Because of the scaling 
properties of TF theory discussed in §2.5, this is not a serious restriction. The results presented 
in §4 are mainly in graphical form but for the convenience of potential users a tabulated version 



can be found at the website http:/ /www.raunvis.hi.is/~ath/TFBT 



Some of the bulk properties of hot dense matter in a strong magnetic field have previously 
been investigated by other authors. Gadiyak et al. (1981) used a modified TF theory including 
more than one Landau level to calculate an equation of state for iron in a magnetic field of order 
10 12 G. Using the strong field version of the TF-method Constantinescu & Moruzzi (1978) studied 
the case of hot iron in bulk in the limit when the field is so strong that only the lowest Landau 
level is populated. A similar investigation was performed by Abrahams &: Shapiro (1991) who also 
constructed a model including exchange and gradient correction to the equation of state of iron in 
the strong field limit. 



2. The Thomas-Fermi Method 

2.1. The General TF Equation 

A quick road to TF theory at arbitrary temperature starts with the basic thermodynamical 
relationship between the particle density n and the pressure P e \ for noninteracting electrons 

n = P' el {n) (1) 

where the dash denotes differentiation with respect to the chemical potential [i. For electrons that 
interact with each other and fixed nuclei by Coulomb forces the electron density depends on the 
position r, and the electrostatic potential generated by K nuclei with charges eZi at positions 
together with the electrons is 

In the TF approximation it is assumed that the relation (|T|) holds locally, but the total 
electrochemical potential /xxf = ^ — e3>n(r) of the electrons (called // t ot in Paper I) is independent 
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of position. The TF equation for the density n(r) thus becomes 

n(r) = F4(e$ n (r) + MTF ). (3) 
The constant //tf depends on the total number of electrons and the boundary conditions. 
The nonlinear integral equation (||) for n may be combined with Poisson's equation 



V 2 $„(r) = 4vre 



A' 



n(r) - - 



i=l 



to obtain the TF differential equation for = + ^tf : 

V 2 *(r) = 4vre 2 P e , 1 (^(r)) 



(4) 



(5) 



for r ^ {R-i, • • • ,Hk}- To take care of the Coulomb singularities at the nuclear positions, (5) has 
to be supplemented by the boundary conditions 



lim |r - Ri|*(r) = e 2 Zi 



r^R 



(6) 



Moreover, the whole system is confined within a bounded region f2, and the total number of 
electrons within $7, N, is fixed. Integrating (||) over $7 we obtain by the divergence theorem 



/ (V*) • dS = Aire 2 ( N - V Z { V 
J an I £W 



(7) 



In the spherically symmetric, neutral case considered in §2.4, 0, is a sphere with radius ro 
containing a single nucleus of charge Z = N at the centre. Hence, the spherically symmetric \E' 
satisfies 



d_ 

dr 



¥(r) 







(8) 



r=ro 



in this case. 



Once (§), (|) and (@) (or (|)) have been solved the electronic density may be obtained from 
(§, i.e. 

n(r) = i^(tf(r)). (9) 
Likewise, the local pressure in the electron gas is 



^ioc(r) = P e l(*(r)). 



(10) 



2.2. TF Theory as a Minimization Problem 

By Legendre transformation one may pass from the pressure P e \ to the free energy density 

/ei(n) = sup{/m - Pel (/•*)}• (11) 
ft 
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The derivative f' el = df e \/dn is the inverse of P' eX , so (||) is equivalent to 

M = fM (12) 

and the TF equation (2) may be written as 

/4(n(r)) - e$ n (r) = Mtf . (13) 

This again states that the total electrochemical potential is independent of the position. Equation 
(|l3| ) is the Euler-Lagrange equation of a minimization problem: Define the free energy functional 
of the interacting electron gas as 

T[n] = f / cl (n(r))d 3 r - £ e 2 Z, f r^-dh + ± [ [ ^^rdV. (14) 

Then is the variational equation for the minimization of under the subsidiary condition 

/n(r )d 3 r = i V = constant . ( 15 ) 



The /iTF in ( J13| ) is the Lagrange multiplier corresponding to the subsidiary condition. It is given 
by 

8F(N) 

where F(N) = ^[utf] is the free energy at the minimizing density nxF> i-e., the solution to the 
TF equation (|i~3|). 



2.3. The Electron Gas in a Magnetic Field at Nonzero Temperature 

We shall now specialize the general scheme described above to the case of non-relativistic 
electrons in a strong, homogeneous magnetic field of strength B. For the rest of this section and in 
§3, we find it convenient to employ atomic units where Planck's constant h, Boltzmann's constant 
k, the electron mass m e and the elementary charge e are all taken to be 1 and dimensionless. This 
means that magnetic fields are measured in the natural unit 

B = = 2.35 x 10 9 G, 

lengths are measured in terms of the Bohr radius 

h 2 o 
a = ^ = 0.529 x 10 -8 cm 

m e e z 

and energies in the unit 

e 2 

— = 27.2 eV. 
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The corresponding temperature unit is 



e 

aok 



and pressure is measured in 



3.16 x 10 5 K 



1 = 8.2 x 10 26 dyn/cm 



(When presenting our results in §4 we shall convert back to conventional units.) 

The motion of the electrons perpendicular to the magnetic field is quantized into Landau 

levels with energy vB, v = 0, 1, 2, The degeneracy of the levels, per unit area, is B/2tt for v = 

but twice as high for v > due to the electron spin. Along the field the motion is one-dimensional 
with the density of states D{e) = e _1 / 2 /(2 1//2 7r), where e is the energy of the translational motion. 

It follows that the particle density at temperature T and chemical potential [i is given by 



n = P e / 1 (/i;T, J B) 



B 1 

27 2V27I- 

BT 1 / 2 



-1/2 



1 



-1/2 



e ( £ +vB- M )/T + ]_ 



2 3 /2^2 



-1/2 



+ 2^I_ 



1/2 



i/=l 



/i — vB 



(16) 



(see figure [l]) where the Fermi-Dirac integral for k > — 1 is 

/o 



4 (x) 



y 

e y-x + i 



Using the relation 4^Ik{%) = klk-i(x) we obtain the pressure 



BT 3/2 



1/2 



i/=l 



(17) 



It is evident that P e i and PL\ depend in a nontrivial way only on the ratios 



V 



T 



and 



B 
T' 



(18) 



Moreover, if £ <C 1 (weak fields and/or high T) and 7? arbitrary, or £ <C |t?| (high or low densities) 
one may use the approximation 



00 roo 1 

£ £ ^(r? - i/O « ^ 4(7? - Ode = ^4+1 



fa)- 



(19) 



Hence, in these limiting cases 



21/2^3/2 



IT" 



1/2 



(20) 
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and 



23/2^5/2 



Pel(»;T)*^r-I 3/ 2{^), (21) 

are independent of B. At low densities, rj — > — oo, one may use the further approximation 



Ik+i(il) ~ (A + l)!e J? , and (20) and (|2l|) reduce to the equation for a classical, ideal gas, 

Pel = nT. (22) 



2.4. TF Theory in the Spherical Approximation 

In the spherical approximation to TF theory of bulk matter each Wigner-Seitz cell is 
approximated by a sphere with the nucleus at the center. Thus, one considers (||), (||) and (0) 
with 17 a sphere of radius ro and a single nucleus with charge Z at Ri = 0. We consider only the 
neutral case, N = Z. Due to the spherical symmetry we may assume that \I/ depends only on the 
radial coordinate r = |r|. To eliminate the nuclear charge Z we define 



x 



Z^r, (23) 



* = T = ^73' (24) 

and it is convenient to replace by 

e W = ^ = f*|^S>. (25) 

(Note that the dimensionless radial coordinate x is not the same as the one used in Paper I 
whereas the dimensionless potential is the same.) 



The TF equation (§) written in terms of these variables becomes 

Pt^x 



d 2 Q(x) _ 2V3 fl _ 1/2 



dx 2 7T 

It has two parameters, (3 and r. 

The approximation (p0|) corresponds to 



7 -l/2 + 2 J -l/2 



XT J ' — ' V XT T 



(26) 



d 2 Q(x) 2 5 /2 , e(a;) 



, o = ' ' ^1/2 

ax z it ' 



XT 



^x h/2 ^± . (27) 



It can be applied when either (3 <C r (high temperature and/or weak fields), or /? <C 0/x (high 
density close to the nucleus). 

The boundary condition @ means that 



6(0) = 1 



(28) 
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and the boundary condition (|8j) is 

e'(xo) = ^1 (29) 

where xq = Z l ^r$, and the dash denotes differentiation with respect to x. 

The pressure Ptf of bulk matter in the spherical TF approximation can be shown to be equal 
to the local electronic pressure at the boundary ( |Fushiki et al. 1989 ): 

^ = ^M = P„,p^) (30) 



and the matter density p is simply given by: 

Am n AZm n 



(47r/3)r-3 (47r/3):cjj 



(31) 



where A is the mass number and m n is the nucleon mass. Combining (|3^) and (|3l| ) gives rise to 
the TF equation of state, 

Ptf = Ptf(p;T,B). (32) 



In practice, one solves ( pq ) with the initial condition ( p8[ ) and a range of values for @'(0). In 
each case one then determines xo from (|29|). By varying O'(0) one obtains the necessary variation 
in the values of xq and hence of p (see e.g. figure [2]). 

The density and pressure profiles in an isothermal surface layer of a neutron star can be 
computed from the equation of state. The surface layer, where TF theory can reasonably be 
assumed to apply in a typical neutron star, is so thin that the gravitational acceleration g s can 
be taken as constant. Denoting by z the depth from the surface, the equation of hydrostatic 
equilibrium for the pressure P is 

dP dPdp(z) 

~a~ = 1 T~ = 9sp ( z >- ( 33 ) 

dz dp dz 

Inserting ( |32|) for P one may obtain p(z) by integration (see a more detailed discussion of 
isothermal atmospheres in §4). As in Fushiki et al. (1989)| and Paper I we find it more convenient, 
however, to express the equilibrium condition in terms of the atomic chemical potential, 

/"atom = FTF + PTFV (34) 

where 

nTF ( r )^3„ , f f n TF (r)n TF (r') 3 3 



Ftf= /el(n TF (rKr - / ^^+1 / ""^"^ ' d*rdV (35) 

Jr<ro Jr<ro T Jr <roJr<ro |r ^ \ 

is the free energy per atom, and 

v = (4vr/3)rg (36) 
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is the volume per atom. The equilibrium condition is flFushiki et al. 1989 ) 

Aiatom^) = ^atom(O) + Am n g s z. (37) 

If Matom is known as a function of the density p, this immediately gives p{z) and hence also P(z) 
from the equation of state. 

In our calculations we have chosen to define the surface of the star to be where the matter 
density, p, matches the lowest density, po, achieved in the zero temperature case. For each case of 
T and B used, we find the value of /iatomO^ = 0) = p a tom(p = Po)- It is then easy to calculate the 
depth z for each value of p from equation ([}?]) and hence evaluate p(z) and P(z). The results for 
isothermal atmospheres are shown in figures 9 and 10 and discussed further in §4. It should be 
kept in mind that free atoms cannot exist at P = if T > 0, so there is no sharp surface. 



2.5. Scaling Relations 

The equations ( |26|) and ( |29"|) have three parameters, (3 = B /Z^l 3 , r = T jZ^I 3 and xq = Z l l 3 r$. 
By scaling they are equivalent to the original equations (|5|) and (||) with four parameters, Z, B, 
T and tq. From (|l7|), ([Til), (p4j) and (|35|) we obtain the corresponding scaling relations for the 
pressure, the free energy and the atomic chemical potential. The matter density p (cf. equation 
(|3lD), is a more convenient variable than ro and we write the scaling relations in the following way: 

P TF (p, B, T; Z', A') = (Z'/Z) 10/3 P TF ((ZA/Z'A')p, (Z'/Z) i/3 B, (Z'/Z) i/3 T; Z, A) (38) 

F TF (p, B, T; Z', A') = (Z'/Z) 10 / 3 F TF ((ZA/Z'A')p, (Z'/Z^B, (Z'/Z^T; Z, A) (39) 

p atom (P, B, T; Z', A') = (Z'/Z) 7 / 3 p atom ((ZA/Z'A')p, {Z'/Zf' 3 B, (Z'/Z) 4 / 3 T; Z, A) (40) 

These relations allow us to compute the pressure, free energy and chemical potential for any pair 
[Z 1 , A 1 ) if it is known for some reference pair (Z, A) as e.g. for iron with Z = 26, A = 56. 



2.6. Limiting cases and validity 

There are three limiting cases, where TF theory simplifies and eventually passes over into the 
theory of a noninteracting electron gas. In terms of the parameters (3 = B/Z 4 ^ 3 , r = T/Z 4 / 3 and 
xq = Z 1 l 3 ro these cases are: 

• High density limit: x (l + /3) 2/5 <C 1. 

• High temperature limit: t/(1 + /3) 2//5 S> 1. 

• Low density limit: rlnxo/(l + /?) 2 ^ 5 » 1. 
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To understand these limits one should recall that for an isolated TF atom at T = the radius tb,z 
is of the order Z" 1 / 3 (l + /3)~ 2 / 5 and the ground state energy E(B, Z) of the order -Z 7 / 3 (l + (3) 2 / b 
( Fushiki et al. 1992| ). The high density limit corresponds therefore to ro <C fB,Zi i-e- highly 



compressed matter. The high temperature limit means that T 3> \E(B,Z)\/Z, i.e. the thermal 
energy is much larger than the ground state energy per electron for isolated atoms. The low density 
limit is mainly of interest as a check for the correctness of the numerical computations. In this 

— 1/3 

limit the electron gas is essentially a classical ideal gas of density n ~ r , because the chemical 
potential ~ Tlnn of free electrons is much lower than the chemical potential ~ E(B,Z)/Z of 
electrons bound to the nuclei. 

In the three limiting cases TF theory may be compared with the simpler uniform model (see, 



e.g., Fushiki et al. 1989 and Paper I), defined by the free energy 



M 9 o / o \ 1/3 

F n (p; B, T) = —U ((Z/M)p; B, T) - -Z 2 (J^j , (41) 

where M = Am n and we use p as a variable instead of rg. The mean electron density is (Z/M)p, 



and the first term in (41) is the free energy of a ball of radius ro of a uniform, noninteracting 
electron gas, whereas the second term is the Coulomb energy of the ball, due to electronic repulsion 
and attraction of the nucleus. The pressure in the uniform model is 

P n (p;B,T) = P el ((Z/M)p;B,T)--(^) Z 2 (J^j , (42) 

with P e i given by (|l7|). Since P e i behaves like p 5 / 3 at high p and as pT at small p or high T it is 



also clear that the first term in (42) dominates over the second in the limiting cases considered. A 



comparison of Ptf with P n and P e \ is shown in figure 6. 



In Hauksson (1996) the status of TF theory with a magnetic field at nonzero temperature 



as an approximation of quantum mechanics is investigated, generalizing the zero temperature 



limit theorems of Lieb et al. (1994) . It is found that thermodynamical functions calculated from 
quantum mechanics converge to the corresponding functions in TF theory in the limit when the 
nuclear charge Z and the magnetic field B tend to infinity, provided B/Z 3 <^ 1. The temperature 
T and the radius ro are in this limit scaled in such a way that r/(l + /3) 2 / 5 and x (l + /3) 2 / 5 arc 
fixed, i.e. T scales as Z 4 / 3 (l + (B/Z 4 / 3 )) 2 / 5 and r as Z _1 /3(l + (P/Z 4 / 3 ))" 2 / 5 . In compressed 
matter and/or at high temperature the condition BjZ 3 <C 1 can be relaxed. In fact by the 
heuristic argument in Paper I, TF theory for compressed matter at zero temperature, i.e. for 
ro < Z~ 1 / 3 (l + (5/Z 4 / 3 )) -2 / 5 , can be expected to be a good approximation, if Z 1 and 
S <C If T 3> B, many Landau levels are excited, and the magnetic field becomes irrelevant. 

Hence at such extreme temperatures the condition for validity is simply Z>1. Summing up, the 
domain of validity of TF theory at non-zero temperature should be at least as large as that for 
zero temperature illustrated in Fig. 1 in Paper I, and strictly larger for T > B. 
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3. Numerical Methods 

Since the calculations presented in this article are an extension of those presented in Paper 
I, the numerical methods are in most aspects similar. The details given in that paper will not be 
repeated, but here we present a general outline together with considerations specific to the case of 
non-zero temperatures. 

To solve the TF equation (|5|) numerically, we use a fourth-order Runge-Kutta method with 
adaptive stepsize control ( Press et al. 1992 ) with the boundary conditions (|28|) and (|29|). First we 



pick a suitable spectrum of values of 0'(O) and then integrate outwards from 0(0) until condition 
( p9|) is satisfied. We then know both xo and Q(xq) which we can use to calculate the density 
and the chemical potential, and hence the corresponding pressure for a given temperature and 
magnetic field. 

It is worth mentioning some changes from the zero temperature case. The kinetic energy 
density, denoted by w in Paper I, is replaced with the free energy density, /. Whereas w is always 
non-negative, the free energy density can become negative because of the added temperature term. 
This allows the electrochemical potential to become negative for high temperatures and low 
densities, since in such cases the free energy, /, can have a negative derivative with respect to n, 
yielding \i < 0. An example of this is seen in figure § where we plot as a function of x for two 
different densities and temperatures with B fixed. Also, when T > there are no free atoms at 
P = 0, contrary to the case for zero temperature (e.g. Paper I). 



As in Paper I, we chose to include the first 100 terms of the infinite sums in the formulas ( |i6[ ) 
and (|l7]) for n and P ensuring a relatively high precision without too much computational cost. 
Here, we encountered a similar divergence problem as in the zero temperature case; close to the 
nucleus, both n and fx tend to infinity. We also found the convergence of the sums in (16) and ( |i~7|) 



to be slow for magnetic fields low compared to the temperature. Both problems may be solved by 
using the approximations ( ^C| ) and (^H), since in the first case £ <C \q\ and in the second £ <C 1. It 
then remains to find a criteria for where to switch over to these approximations. It is convenient 
to rewrite the expressions depending on [i in equations (|l^) and (|i~7|) as 



4(£C) + 2]T4(£(C-^) 



u=l 



with 



S B 



and determine in terms of £ and £ where to switch over to the approximations to obtain the desired 
accuracy. We decided that in the worst case the relative error in P and n should never exceed 
2 x 10~ 4 for any value of fi at a given temperature and magnetic field strength. We calculated 
this error for the two methods, using only 100 Landau levels on one hand and the approximations 
on the other and comparing the results with the values obtained with 5000 Landau levels (which 
we consider as exact). This comparison yielded critical values for f, referred to as £ c , above which 
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the approximations ( poj) and pl| ) are used instead of the partial sum over 100 levels (see table [j]). 
Since different Fermi-Dirac integrals enter the calculations of n and P, slightly different £ c values 
are appropriate in each case for the same £. The Fermi-Dirac integrals were evaluated using the 
method of Cody & Thatcher (1967) which assures a relative error less than 3 x 10~ 9 . (Note that 
in their paper, the value of q s for s = and n = 4 in table IIC should read 1.00 and not 1.05.) 



4. Results 

In this section we present the results of our numerical TF calculations for bulk matter in high 
magnetic fields and at finite temperature. We remind the reader that our reference element is iron 
with Z = 26 and A = 56 and that the treatment is non-relativistic. Hence the highest density 
used in our calculations is 10 6 g cm -3 and the highest temperature is 100 keV (1.16 x 10 9 K). 
As in Paper I the range in magnetic field strength is from 10 10 to 10 13 G. The reason for this is 
that in the case of iron one obtains the zero-field results (see e.g. Abrahams fc Shapiro 199l| ) if 



the fields are much weaker than 10 10 G whereas for field strengths slightly higher than 10 12 G 



essentially only the lowest Landau level is populated and one obtains the results of Constantinescu 



fc Moruzzi (1978)| and Abrahams fc Shapiro (1991) 



Our presentation of results is mainly in graphical form but tabulated results are available 



at the website http: / / www.raunvis.hi.is / ~ath /TFBT . By use of the scaling relations discussed 



in §2.5 the numerical TF results can also be used for elements other than iron and for different 
temperatures and magnetic field strengths. 

Figure |3| shows a typical example of how temperature changes the electron distribution inside 
a unit iron cell of a given size, corresponding to fixed matter density The quantity plotted is 
r 2 n(r) which is useful since the the number of electrons in a shell is proportional to the area 
under the corresponding curve. The case shown is for a field strength of 10 11 G and a density of 
1000 g cm -3 . The main thing to notice is that thermal effects tend to smooth out the distinct 



Landau shell structure which is so prominent in the zero temperature case (see Fushiki et al. 



1992, and Paper I). However, as long as the temperature is lower than approximately one tenth 
of the cyclotron energy, Tiuj B (« 11.58Si2 keV where B 12 = B/W 12 G), thermal effects are not 
particularly pronounced. It is only when the temperature is comparable to or higher than the 
cyclotron energy (~ 1 keV for 10 11 G) that the electron distribution becomes drastically different 
from the zero temperature case. This is a reflection of the temperature dependence of the relation 
between n and [i for the free electron gas shown in figure |l]. From figure |3| it is also clear that the 
number density of electrons at the edge of the unit cell increases with increasing temperature and 
thus the TF pressure of bulk matter at given matter density and field strength also increases with 
temperature. 

In figures 4 to 10 we present the results of our calculations of the bulk properties of dense 
matter. Figure 4 shows the TF equation of state for iron, Ptf versus p, where Ptf is given by 
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equation (^) and p by equation (^T]). Each curve corresponds to a fixed value of the temperature. 
The presence of a Landau shell structure in the electron distribution is reflected in a clear 
oscillatory nature of the pressure-density relationship as long as the temperature is lower than 
approximately one tenth of the cyclotron energy (i.e. T < l-B^keV). For somewhat higher 
temperatures the oscillations are smoothed out as previously noted by |Gadiyak et al. (1981) . 



In figure ^c] we compare the TF equation of state for a field strength of 10 12 G with that of a 
nondegenerate gas of free electrons with uniform density n = Zp/Am n and temperature T but in 
zero magnetic field. This pressure is given by equation (|22] ) or equivalently by 



^Pciass = PdueiP, T,B = 0) = ZpkT/Am n , (43) 

where we use the notation P c iass instead of P e \ m order to avoid confusion with the more general 
expressions for P e \ in §2.3 (also note the change of units from those used in §2). We have chosen 
the value of 10 12 G for this particular discussion but the same qualitative features are present for 
all B. Figure ^ shows that at low densities the magnetic TF pressure approaches -Pciass for fixed 
B and T. The same is true in the limit of high temperatures at fixed p and B. Hence at very low 
densities and/or very high temperatures the TF pressure becomes independent of B as expected 
from the discussion of asymptotic behaviour in §2.6. On the other hand at high enough densities 
and as long as the temperature does not become too high both thermal and magnetic effects on 
the equation of state become negligible and it approaches the non-magnetic TF equation of state 
of iron at zero temperature (Paper I). This limiting behaviour is clearly seen in figure || which 
shows Ptf versus p for T = 0.1 keV and several values of the magnetic field. 

In figure 6 we compare the TF pressure, Ptf, with P n , the pressure of the uniform model 
discussed in §2.6 (see equation (fi"2"|)). Results are shown for a magnetic field with strength 10 12 
G and two values of the temperature. Also shown is the pressure of a noninteracting electron 
gas, P e i, given by equation (|17]). It is clear from the figures that the three models have the same 
behaviour in the limit of high or low densities and at high temperatures as already discussed in 
§2.6. The differences between the different approximations are greatest at low temperatures and in 
a restricted range of relatively low densities. This is approximately the range in which one expects 
contributions from the thermal motion of the nuclei to be important (see the discussion below). 

For completeness we show in figure 7 the Thomas- Fermi free energy for iron in bulk, -Ftf> as 
a function of p for selected values of T and B (see equation (||) ) . 

At this point we would like to mention a few corrections to the equation of state that go 
beyond the TF approximation and are therefore not included in our calculations. For further 



discussion of this and related topics we refer the reader to Paper I and the papers by |Fushiki et al. 
19921 and|Yakovlev fc Kaminker (1994) 



As already mentioned in the introduction the pure TF method neglects both exchange effects 
between the electrons as well as "gradient corrections" of the von Weizsacker type. Although 
these effects are not expected to make much difference as far as bulk properties of matter are 
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concerned ( [Fushiki et al. 1992 ) it should be kept in mind that they may become important at very 



low temperatures and the densities corresponding to the surface densities of neutron stars ( |Fusliiki 



et al. 1989; Abrahams fc Shapiro 1991] ; Paper I). To investigate this in detail requires more exact 



calculations than are provided by the one dimensional (spherical) TF method. 

An additional simplification is that in TF theory the nuclei are assumed to be stationary 
at fixed relative distances determined by the matter density, their number density given by 
n; = p/Am n . Thus the thermal motion of the nuclei is neglected in calculations of the TF pressure. 
This motion is only weakly affected by the magnetic field and for most applications of relevance to 
the investigation of neutron star envelopes it is a fairly good approximation to regard the nuclei 
as forming an ideal non-degenerate gas with pressure 

Pi = mkT = pkT/Am n (44) 



(for a more detailed discussion see e.g. Gudmundsson et al. 1983| , and Yakovlev & Kaminkci 



|1994 ). Note however that Pi = P c \as S /Z, where -P c iass is the pressure of the ideal non-degenerate 
gas of electrons at the same p and T (see equation (|i"3|)) and thus the pressure of the nuclei is 
small compared to the TF pressure except at low temperatures and and low densities, which are 
of little importance for compressed matter. This can e.g. be seen by inspecting figures ^ and 6 
which show a comparison between PrF) Pciassj Pel, an d P u for B = 10 12 G. We thus feel justified in 
neglecting the thermal motion of the nuclei in our calculations. We point out however that in case 
one wants to include the contribution of the nuclei to the total pressure and chemical potential 
this can easily be done by adding them to the TF results, thus obtaining 

P = P TF + P; (45) 

and 

V = ^atom + W , (46) 

where /U at0 m is given by equation (^4|) . For an ideal non-degenerate gas of nuclei in zero magnetic 
field Pi is given by equation ([|4|) above and 

Pi = kT In [ni(2-Kh 2 /m n kT) 3/2 ]. (47) 

A more detailed treatment taking into account the non-ideal behaviour of nuclei due to Coulomb 
interactions between nuclei and quantum effects is beyond the scope of this paper. We would 
however like to remind the reader that these effects are expected to make only relatively small 
corrections to the equation of state of bulk matter for most of the temperatures and densities of 
interest in the study of neutron star envelopes. On the other hand they do have a considerable 
influence on the transport properties in magnetic neutron stars, and they also determine which 



form or phases of matter are present in the envelope (see e.g. Gudmundsson et al. 1983 ; Van Riper 



1988 : Yakovlev fc Kaminker 1994 ; Lai Sz Salpeter 1997] , and references therein). 
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We now return to the discussion of the results of our TF calculations. Figure |8| shows the 
isothermal sound velocity, 



a T = — (48) 



dp ) 



T 



with P = PtFj as a function of density for several values of the temperature and a magnetic field 
of 10 12 G. This quantity is related to the adiabatic sound velocity a s = (dP/dp)^ 2 by 



2 c v 2 

a T = —a s , 



where c p and c„ are the specific heats at constant pressure and volume, respectively. As discussed 
in Paper I (see also Landau fc Lifshitz 1960| ) a s is the velocity of longitudinal magnetosonic waves 



when the wavevector is parallel to the magnetic field. Figure [8| shows that the oscillatory behaviour 
which is due to the Landau shell structure in the electronic density distribution of the unit cells 
is clearly present in aT at finite temperatures. Even though thermal effects tend to smooth them 
out they are quite pronounced as long as the temperature is less than approximately one tenth 
of the cyclotron energy. It is only when the temperature is of order 10 B\2 keV or higher that 
the oscillations disappear completely, an effect that has already been discussed in connection with 
the equation of state. In two or three dimensional TF calculations an additional smoothing effect 
would probably be present due to the non-spherical nature of the unit cells. However as argued in 
Paper I they would not be big enough to wipe out the oscillations. 

These de Haas-van Alphen type oscillations are reflected in the isothermal compressibility of 
bulk matter Kbuik, which is related to the sound velocity by 

Kbulk = —^o- ! ( 49 ) 
pa T 

and hence they are also present in the density and pressure profiles of matter in the surface layers 
of neutron stars (Paper I). In most cases of interest the temperature in the surface layers rises 
quite steeply as a function of depth and hence the oscillatory behaviour in the matter density can 
not be investigated in detail without solving the thermal structure equations of the envelope. This 
requires knowledge of the opacity as well as the equation of state ( pudmundsson et al. 1983| Van 



Riper 1988| ). Some insight may clearly be gained however by considering the simplified case of 



isothermal surface layers, an approximation which may in fact be appropriate for sufficiently old 
and relatively cold neutron stars. Figures 9 and 10 show the results of such isothermal structure 
calculations for a range of temperatures and two values of the magnetic field strength. The method 
used in the calculations is that described at the end of §2.4. The oscillatory behaviour is clearly 
seen in figure 9 as steps in the density profiles for the lowest temperatures whereas the effect is 
somewhat weaker in the pressure profiles in figure 10. Notice also that at high temperature both 
the density and pressure gradients are very small in the outermost layers. This is a result of the 
fact that there is no sharp boundary at T > 0. For this reason the depth shown in figures 9 and 10 
is measured from the point where the surface would be at zero temperature (see the discussion at 
the end of §2.4 and Paper I). 
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In some applications it is also of interest to know the column density, a, either as a function 
of zg s or of p. In either case it is easy to calculate the column density from the surface to the point 
in the star where the pressure is P since it is simply given by the expression 

f z f p dP' 

a-a s = pdz' = = (P - P s )/g s , (50) 

Jo Jp s 9s 

where we have used equation (^3[) and a s and P s are the column density and pressure at the 
surface of the neutron star, respectively. Hence the results shown in figures 4 and 10 can easily be 
used for evaluating a. 

We end the discussion of our results by estimating the effects of the thermal motion of the 
nuclei on the velocity of sound and hence also on Kbulk ana - the isothermal surface structure. If we 
assume as before that the nuclei form an ideal non-degenerate gas at temperature T then by use 



of equations (44) and (E3) we find that the TF sound velocity shown in figure and denoted by 



ar has to be replaced by a T given by 

a' T = (a 2 T + kT/Am n ) 1/2 . (51) 

By comparing the nuclear contribution, {kT / Am n ) 1 ^ 2 ~ 4 x 10 6 (T/keV) 1//2 cm/s, to the values 
of Or shown in figure || one can see as before that it is only at the lowest temperatures and in 
a restricted range of low densities that the thermal motion of the nuclei needs to be taken into 
account. In particular it is clear that this motion has a negligible effect on the high density 
oscillatory behaviour which occurs at low temperatures. 



5. Conclusions 

In this paper we have presented the results of a detailed calculation of the properties of hot 
dense matter in high magnetic fields using the temperature dependent Thomas-Fermi method and 
taking all Landau levels into account. By including temperature these calculations extend our 
previous work on the properties of cold dense matter in the TF approximation presented in Paper 
I. 



On the grounds discussed in [Fushiki et al. (1992)| and Paper I we expect that TF theory 



catches the main features of bulk matter composed of heavy elements like iron. This applies also 
for finite temperatures and high magnetic fields up to at least a few times 10 12 G. 

We have investigated the range of validity of the finite temperature TF model for matter in a 
magnetic field and estimated the effects of the thermal motion of nuclei on the properties of matter 
in bulk. We find that the nuclear contribution to the total pressure and other bulk properties is 
very small except in a restricted range of low densities and low temperatures. 



Within the range of validity of the finite temperature TF theory we have found that the 
Landau shell structure which is so pronounced in iron at zero temperature is also present at 
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finite temperature as long as it does not exceed about one tenth of the cyclotron energy. This 
structure is reflected in an oscillating component of the de Haas-van Alphen type in all the 
bulk properties of matter as functions of depth in isothermal surface layers of magnetic neutron 
stars. At temperatures around the cyclotron energy or higher, thermal effects smooth out the 
oscillations. Further investigation of this problem including the additional effects of smoothing 
of oscillations due to the non-spherical shape of the unit cells requires two or three dimensional 
TF calculations. Furthermore an investigation of the thermal structure of the surface layers when 
the assumption of isothermality is not valid requires solving the thermal structure equations of 
neutron star envelopes. This in turn requires knowledge of the transport properties of matter 
in high magnetic fields in addition to the equation of state presented in this work. Such study 
is beyond the scope of this paper. Finally we remind the reader that our numerical results are 



available in tabular form at the website http:/ /www.raunvis.hi.is/~ath/TFBT 
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B/T 


0.0001 


0.001 


0.01 


0.1 


1 


10 


100 


n 






32 


97 


97 


97 


97 


P 




50 


97 


97 


97 


97 


97 



Table 1: Values of £ c . For £ > Cc approximations ( |20| ) and ( |2~l| ) are used when calculating n and 
P, respectively. Where no value is given, the approximations are used for all values of (. Here, B 
is measured in units of 10 12 G and T in keV. 
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Fig. 1. — The electron number density as a function of the chemical potential for a uniform gas 
of free electrons. The curves are drawn for B = 10 11 G and three temperatures between 0.01 keV 
and 1 keV in addition to the zero temperature case. Note that the Landau band structure which 
is still present at T = 0.01 keV has almost disappeared at T = 0.1 keV and is not visible at T = 1 
keV. The quantities on the axes are dimensionless: n is written in units of a~ 3 and fi in units of 
the cyclotron energy (equal to B in dimensionless units). 
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Fig. 2. — The potential as a function of the radial coordinate x for two unit cells corresponding 
to two different values of the matter density. Results are shown for two values of temperature and 
a magnetic field strength of 10 12 G. The behaviour of for T = 0.01 keV is very similar to the zero 
temperature case. Note that for higher temperatures, such as T = 1 keV shown here, 6 becomes 
negative at the low matter density whereas at the higher density it takes values which are higher 
than in the zero temperature case. 
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Fig. 3. — The quantity r 2 n for a spherical unit cell of iron in the TF approximation. The matter 
density is 1000 g cm" 3 and the magnetic field strength 10 11 G. Results are shown for five different 
values of the temperature. 



November 6 1997 



24 




p [g/ cm3 ] 

Fig. 4a. — The Thomas-Fermi equation of state for iron in a magnetic field. In each case results 
are shown for six values of the temperature, including the zero temperature case, (a) B = 10 10 G. 
(b) B = 10 11 G. (c) B = 10 12 G. (d) B = 10 13 G. In 4c the pressure of an ideal non-degenerate 
gas of free electrons in zero field, -P c iass, is superimposed on the TF results for comparison. See the 
text for further explanation. 
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Fig. 4b.— 
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Fig. 4d.— 
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Fig. 5. — The Thomas- Fermi equation of state for iron at a temperature of 0.1 keV and several 
values of the magnetic field strength. Note in particular that the curve for B = 10 10 G falls almost 
completely on top of the zero field curve. 
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icr 10' icr icr 10" icr icr 

p fg/cm 3 l 

Fig. 6a. — The equation of state for iron in three different approximations. The TF pressure, -Ptf, 
is shown by the solid line. The pressure of a Fermi gas of free electrons, P c \, is shown by the dashed 
line and the pressure in the uniform model, P u , by the dotted line. The strength of the magnetic 
field is 10 12 G. (a) T = 0.1 keV. (b) T = 1 keV. 
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Fig. 7a. — The free energy of iron in bulk, Ftf, in the Thomas- Fermi approximation. Results are 
shown for several values of the temperature, (a) B = 10 10 G. (b) B = 10 12 G. 
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Fig. 8. — The isothermal sound velocity in bulk matter made of iron in the Thomas-Fermi 
approximation. Each curve corresponds to a given value of the temperature as indicated. The 
inset shows details around the peak of the first oscillation. 
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Fig. 9a. — The density of isothermal surface layers of neutron stars as a function of depth in 
the Thomas-Fermi approximation. Results are shown for five values of the temperature and two 
values of the magnetic field strength. The surface gravity is written in units of 10 14 cm s~ 2 . Since 
isothermal atmospheres extend to infinity in the TF approximation the depth shown is measured 
from the point where the surface would be at zero temperature, (a) B = 10 10 G. (b) B = 10 12 G. 
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Fig. 9b.— 
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Fig. 10a. — The pressure of isothermal surface layers of neutron stars as a function of depth in the 
Thomas-Fermi approximation. Results are shown for five values of the temperature and two values 
of the magnetic field strength. See the text and the caption for figure 9a for further explanations, 
(a) B = 10 10 G. (b) B = 10 12 G. 
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Fig. 10b.— 



